Geodesies and the best measurement for distinguishing quantum states 
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From statistical distinguishability of probability distributions, one can define distinguishability 
of quantum states. A corresponding measurement to perform, optimal in a definite sense, for 
distinguishing between two given states pA and ps, has been derived by Fuchs and Caves. We show 
that the Bures-Uhlmann geodesic through the two states singles out this measurement. The geodesic 
'bounces' at the boundary of the set of quantum states. Whenever the geodesic hits the boundary, 
the state orthogonal to that boundary state is one of the basis states for the measurement. 



Encountering two quantum states one might ask how 
similar they are. Clearly one can think of ways to un- 
derstand 'similar' that would make this question relevant 
in quantum information processing. One reformulation 
of the question could be: How well can we distinguish 
between the two states, with the aid of a measurement? 
Orthogonal states are one-shot distinguishable — it is pos- 
sible to measure just once to know the state. But for 
given non-orthogonal states there is no measurement that 
will discriminate between the states with certainty. Then 
we can ask for the measurement that will be the most 
advantageous one, in a statistical sense, for distinguish- 
ing them. For every measurement, the two states give 
two probability distributions for the outcomes. To distin- 
guish between the states we need to distinguish between 
the probability distributions. 

This leads to statistical distance [II P| between states 
and, in a definite sense, a best distinguishing measure- 
ment. Fuchs and Caves [3 have derived an expression for 
this distance, and also for the operator corresponding to 
the optimal measurement. The distance turned out to be 
the same as the geodesic distance of the Bures-Uhlmann 
metric 0,0,00 set of quantum states. In this let- 

ter we first introduce these results and some facts about 
the Bures-Uhlmann metric. This will lead the way to our 
result. We will show how the optimal measurement is de- 
termined by the Bures-Uhlmann geodesic connecting the 
two quantum states: the geodesic 'bounces' at the bound- 
ary of the set of quantum states, at states orthogonal to 
the basis states of the measurement operator. 

Every measurement that can be performed on a quan- 
tum system can be described by a POVM — a positive 
operator valued measure. This is a complete set of 
non-negative Hermitian operators Ei acting on an A'^- 
dimensional Hilbert space TLn] * indexes the measure- 
ment outcomes and completeness means Ei = 1 (the 
identity operator). If the POVM-elements are orthogonal 
one-dimensional projectors, Ei = |ei)(ei|, we have a von 
Neumann measurement, with a corresponding observable 
O — ^f^^ XiEi. Upon measuring, the outcome i will oc- 
cur with probability pi = TrEip, where p is the density 
operator describing the state of the system. Thus, two 
states Pa and pB will produce two probability distribu- 
tions p*^"*^ and p(^^ 



Sampling from these distributions a finite number of 
times will give frequencies that differ somewhat from the 
probabilities. From the size of the statistical fluctua- 
tions, a measure of distinguishability between probabil- 
ity distributions can be found. This is the statistical 
distance d{p^'^\p^^^), by Bhattacharyya and Woot- 
ters 0, given by 

co.dip(^\pi^y)^j2^M^^ ■ w 

i 

The corresponding Riemannian metric is known as the 
Fisher-Rao metric. 

Different POVMs give different probability distribu- 
tions, and thus different statistical distance. Take the 
measurement that gives the maximal statistical dis- 
tance. Then, by definition, that distance is the distance 
d{pA,pB) between quantum states: 

d{pA:PB) = max arccos ^ ^jTrE^^jTrE^p^ . 

_ (2) 

Fuchs and Caves [3 showed that the maximization yields 



d{pA,PB) = arccos Try p^^ps pj/^ , (3) 

where ^ = p^/^ is the unique non-negative operator 
such that p^/'^p^/'^ = p. This expression looks asymmet- 
ric, but one can show that it is not. The most separated 
states are orthogonal states, which rest at the distance 
tt/2 from each other. For pure states \iPa) and \iPb) we 
have d — arccos the angle in Hilbert space 

between the two vectors. 

It may be noted that the same trace expression appears 
in the fidelity 0: 

F{pa,Pb) ^(rr^ pTpb pT) ■ (4) 

It is a choice for a quantity that should measure the ac- 
curacy of transmission of a signal via a quantum channel. 

What is more of our concern here is that the distance 
d{pA,PB) ijSJ is the same as the geodesic distance in the 
set of density operators, according to the Bures-Uhlmann 
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metric. Shortly we will see how this metric is constructed. 
But first we will take a look at the best measurement — 
the measurement found to correspond to the maximal 
distance (|3J) 0] . It is a measurement of the observable 



1/2 



1/2 1/2 -1/2 

Pa Pb Pa Pa 



(5) 



This optimal measurement is unique, except for some 
special cases. (These are when M has degenerate eigen- 
values, or when the states are non-invertible, i.e. when 
they are boundary states.) Although we have an explicit 
expression for the best distinguishing measurement, it is 
not easy to determine what it is for given density opera- 
tors. Operator square roots are rather unwieldy to com- 
pute. However, this operator has appeared earlier. It is 
a part of the prescription for Bures-Uhlmann geodesies. 
We will take advantage of this connection to get a new 
description of the best measurement. 

It should also be mentioned that formula jSJ is an oper- 
ator mean. Operator means, fulfilling a set of reasonable 
criteria, can be defined for two positive operators A and 
B 0. The geometric mean. 



Ai^B = A^/\ A-y^BA-y^ 



(6) 



is one of these. Thus, we see that the optimal observ- 
able lO is the geometric mean of p^^ and pB'- 



(7) 



M = p/#pB 



Now we will turn to the Bures-Uhlmann metric 
0, 0|. It is obtained from a (kind of) fibre bundle con- 
struction in the Hilbert-Schmidt space. From this space 
of operators VF, acting on Tijv, we have the projection 
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WW'' 



(8) 



to the base manifold of positive operators. {W'' denotes 
the Hermitian conjugate of W.) The fibres are obtained 
by right multiplication of the unitary group, since W and 
WU, for unitary U, will be projected to the same opera- 
tor. In the bundle space we define distances D{Wa, Wb) 
by 

D^{Wa,Wb) = Tr(W^^ - Wb){w\ - wj,) (9) 

-a Euclidean distance. In the base manifold of positive 
operators, the distance between two operators is defined 
as the length of the shortest path between the corre- 
sponding fibres in the bundle space. Here we are only in- 
terested in the set of normalized density operators, which 
means projections of operators W on the unit sphere, 
TrVFVFt = 1, in the Hilbert-Schmidt space. With this 
restriction, the geodesic distance we get in the set of den- 
sity operators is given by equation (j^J. 

The physical interpretation of this construction is that 
of state-purification. Every mixed state p can be puri- 
fied in a larger Hilbert space; the system is regarded as 



s subsystem of a bipartite system. The Hilbert-Schmidt 
space takes the role of the larger state space. Every op- 
erator W represents a pure state vector, and the reduced 
density operator for the subsystem is WW^ . W is said 
to be a purification of p, if p = WW^ . Consequently 
the whole fibre WU consists of purifications of p. The 
distance between the states pA and pb of the subsystem 
should not be larger than between any two purifications 
Wa and Wb- This is assured by finding the shortest path 
between the fibres of purifications . 

Geodesies on the unit sphere (i.e. great circles) in the 
Hilbert-Schmidt space can be expressed in the following 
way; 

W{t) = WoCOst + Woaiiit , < i < 27r , (10) 

with normalization and orthogonality conditions, 
TtWoW^ = TiWoW^ = 1 and Tr(W^oW^cf + W^Wq) = 0. 
For this curve to project to a geodesic in the set of 
density operators, it is also required that it is every- 
where perpendicular to the fibres — the 'horizontality 
condition'. This reads 



Wr'Wn = W'W< 



^^0 



(11) 



The geodesic of density operators is then the curve 

p(t) ^W{t)W''{t) , 0<t<TT. (12) 

The projected curve will do two turns, since W 
and —W — always resting on the same geodesic — will be 
projected to the same density operator. Hence, the range 
for t is halved. 

Let us consider the geodesic between the states pA 
and Pb- Assume Wa to be a preimage of pA- Using 
the 'horizontality condition' Hll|l (and positivity of pa), 
it can be shown that the operator Wb, that is, the preim- 
age of Pb, should be given by 



Wb = MWa , 



(13) 



where M is the positive operator given by jS} — the op- 
erator corresponding to the best distinguishing measure- 
ment. 

Now we start the geodesic at pA (assumed to be in- 
vertible) and let it go through pB, resting at the geodesic 
distance d away. This means that we set 

(Wa = WiO) , 

I Wb — W{d) , where cosd = Tr^J pIC^Pb Pa^ ■ 

(14) 

Inserting in (|10|l . we can solve for the geodesic W{t) in 
terms of Wa and Wb ■ 

sint 



W{t) = Wa cost -I- {Wb - Wa cosd) 



sind 



(15) 



Alternatively, if we use Wb — MWa, we can express 
W{t) in terms of Wa and M: 



Wit) = lcost+{M - Icosd)- 



s'mt 



sind 



Wa 



(16) 
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Here d should be understood as given by cosd = TrMpA- 
Finally we have, for the projected curve [lO| . 

p{t)=X{t)pAX{t) , 

where X{t) = (Icost + {M - ]lcosd)f^) . 

In this form we can think of a geodesic as given by a 
starting point pA and a positive matrix M, determining 
the direction from pA- 

This shows that there is a close relation between the 
geodesic and the operator M for the optimal distinguish- 
ing measurement. It is then natural to ask: To what ex- 
tent does the geodesic determine the best measurement? 
First we note that the measurement is given by the eigen- 
basis of M, while the eigenvalues are superfluous infor- 
mation. To answer the question we need to know more 
about the geodesies. In one of Uhlmann's papers 01 it 
is explained how the Bures-Uhlmann geodesies 'bounces' 
at the boundary of the set of density operators. We will 
now investigate this feature with the purpose of proving 
a geometric description of the optimal measurement. 

At the boundary of the set of quantum states the ma- 
trices p have at least one zero eigenvalue, hence det p = 0. 
Consider the states on a geodesic, given by equation (fT7|l . 

det p{t) = 4^ det = 

^ det {M - Ix) = , (18) 
where x — cos d — cot t sin d 

The solutions for x of this equation are the eigenvalues 
Aj of M. For the corresponding values of t, the states 
p{t) lie on the boundary. These states are 

sin"^ f ■ 

p{t,)^{M~lX,)pA{M-lX,)-^ . (19) 

sm d 

{ti is given by the equation — cos d— cot U sin d.) There 
are iV or less boundary states, since it is the same as the 
number of different eigenvalues of M . And the number of 
zero eigenvalues of p(ti) is the same as the degeneracy of 
the relevant eigenvalue of M. When the projected curve 
reaches the boundary of the set of density matrices it will 
bounce back into the interior. After N, or sometimes less, 
bounces, the curve will return to its starting point. Thus, 
the geodesic will consist of a number of segments with 
endpoints at the boundary of the set of density matrices. 

From the expression H19|l . it is easily recognized 
what states are orthogonal to geodesies boundary 
states p{ti) — they are nothing but the eigenvectors \mi) 
of M: 

{m,\p{U)\m,) = . (20) 

If all eigenvalues A^ are distinct, we get N boundary 
states, each with one zero eigenvalue, which singles out 
the N basis states \mi). If there are degeneracies we get 



less boundary states, but the sum of the zero eigenval- 
ues is stiU N. And a state p{ti) with n zero eigenvalues 
is orthogonal to an n-dimensional subspace of the pure 
states. In this subspace any basis can be chosen; M is 
diagonal in anyone of them. Thus, the geodesic singles 
out the basis states — the optimal measurement is fully 
determined by the geodesic. 

For two-level systems this provides a practical method 
for finding the optimal measurement. Every pair of states 
then lies in a disc in the Bloch ball, which is isometric to a 
round hemisphere \vi\ . The geodesies on the hemisphere 
are just great circles and the relation to the disk is by 
orthographic projection. In this case the geodesic's end- 
points at the 'equator' are the basis states of the measure- 
ment. For higher dimensional Hilbert spaces, the picture 
is much more complex. 

In conclusion we have seen how the Bures-Uhlmann 
geodesies in the set of quantum states bounces at the 
boundary in a set of N, or less, states. If there are N 
boundary states, the states orthogonal to these form an 
orthonormal basis in the Hilbert space. If there are less 
than boundary states there are just enough lower rank 
states, so that it is again possible to form an orthonor- 
mal basis of states orthogonal to the boundary states. We 
have shown that such a basis is the basis of an optimal 
distinguishing measurement: the best measurement — in 
a specific statistical sense — to perform for distinguishing 
between two quantum states lying on a segment of the 
geodesic. This result is a conceptually interesting char- 
acterization of the optimal distinguishing measurement. 

I gratefully acknowledge many discussions with Inge- 
mar Bengtsson. He has taught me about the Bures- 
Uhlmann geometry and encouraged me to do this work. 

Note added in proof: An expression for the geodesies, 
similar to equation (|17(l . has also been derived by Bar- 
num [13. 
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